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Abstract 


A serpentine flow channel is one of the most common and practical channel layouts for a polymer electrolyte membrane (PEM) fuel cell since 
it ensures the removal of water produced in a cell with acceptable parasitic load. During the reactant flows along the flow channel, it can also leak 
or cross to neighboring channel via the porous gas diffusion layer due to the high pressure gradient caused by the short distance. Such a cross flow 
leads to a larger effective flow area altering reactant flow in the flow channel so that the resultant pressure and flow distributions are substantially 
different from that without considering cross flow, even though this cross flow has largely been ignored in previous studies. In this work, a numerical 
and experimental study has been carried out to investigate the cross flow in a PEM fuel cell. Experimental measurements revealed that the pressure 
drop in a PEM fuel cell is significantly lower than that without cross flow. Three-dimensional numerical simulation has been performed for wide 
ranges of flow rate, permeability and thickness of gas diffusion layer to analyze the effects of those parameters on the resultant cross flow and 
the pressure drop of the reactant streams. Considerable amount of cross flow through gas diffusion layer has been found in flow simulation and 
its effect on pressure drop becomes more significant as the permeability and the thickness of gas diffusion layer are increased. The effects of this 
phenomenon are also crucial for effective water removal from the porous electrode structure and for estimating pumping energy requirement in a 
PEM fuel cell, it cannot be neglected for the analysis, simulation, design, operation and performance optimization of practical PEM fuel cells. 


© 2006 Elsevier B.V. All rights reserved. 
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1. Introduction 


Use of fossil fuels and consequent damage to the atmosphere 
has become one of the most serious concerns of modern times, 
and has called for environmentally friendly energy conversion 
and power generation systems for sustainable development. 
Polymer electrolyte membrane (PEM) fuel cells are a clean 
power source, at least at the point of use, and have been regarded 
as one of the most promising advanced energy technologies. 
Significant effort is being spent worldwide to develop PEM fuel 
cells for mobile, portable and stationary applications. 

The performance of a PEM fuel cell is influenced immensely 
by the flow distributions in the flow channels and gas diffusion 
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layer (GDL) since water flooding and concentration distribution 
are the most important issues governed by the configurations of 
flow channel layouts. Watkins et al. [1] proposed a continuous 
fluid flow channel that had an inlet at one end and an outlet at the 
other, and typically followed a serpentine path. Such a serpen- 
tine flow channel layout can be regarded as many parallel flow 
channels being connected in series. This flow channel design 
results in a long reactant flow path; consequently, a significant 
pressure drop occurs with considerable concentration gradients 
from the channel inlet to the outlet. Hence several separate flow 
channels are often adopted for large cells in order to reduce the 
pressure drop and excessive power requirement for the supply 
of air stream [2]. Classifications and details of discussions on 
flow field designs can be found in Li [3] and Li and Sabir [4]. 
Serpentine flow channel layout and many of its modifications 
are widely chosen as flow channel designs for PEM fuel cells 
owing to the fact that this flow channel design ensures effec- 
tive water removal from the cell, while with reasonable pressure 
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Nomenclature 


Acell active cell area 


Cr wall friction coefficient 

C a constant in the expression for the laminar wall 
friction coefficient 

d flow channel depth (m) 

dh flow channel hydraulic diameter (m) 

F Faraday constant 

H bipolar plate effective height (m) 

J cell current density (A cm~?) 

Len flow channel length (m) 

Neurn number of turns in flow channel 

P pressure 

Q volume flow rate (m? s7!) 

R universal gas constant (8.314kJ kmol~! K7!) 

Ro outer diameter for U-turn of flow channel 

Re reynolds number 

S distance between flow channels (land width) 
(mm) 


T temperature (K) 

u velocity vector (m s7!) 

V average velocity (m s7!) 

w flow channel width (mm) 

W bipolar plate effective width (m) 


Greek letters 


A difference 

ôe thickness of gas diffusion layer 
H viscosity (N s m~?) 

p density (kg m~?) 

¢ stoichiometry 

Subscripts 

cell fuel cell 

e electrode (gas diffusion layer) 
f friction 

GDL gas diffusion layer 


inlet value at the inlet of the flow channel 
max maximum value 


min minimum value 
outlet value at the outlet of the flow channel 
turn U-turn 


drop, resulting in excellent cell performances. The gas diffusion 
layer (GDL) also has considerable effects on the flow distribution 
in a PEM fuel cell. The thickness, porosity, permeability and the 
wetting characteristics of the pores are most important parame- 
ters which define the characteristics of the GDL. The pores can 
be hydrophobic with optimal value of polytetrafluoreothylene 
(PTFE) [5] not to be congested with liquid water. Water flow 
in GDL is experimentally investigated for various thicknesses 
and pore sizes by Benziger et al. [6]. Paganin et al. [5] observed 
a performance decrement at higher current densities when the 
diffusion layer thickness is increased. The thickness of the GDL 
was optimized by a mathematical modeling in Inoue et al. [7] 


and by cell performance tests in Lee et al. [8] and Jordan et al. 
[9], with various GDL parameters such as porosity, permeability 
and thickness. A thin GDL with small porosity results in good 
electrical conductivity, however efficient mass transport requires 
large pores. 

Multi-dimensional numerical simulations have been carried 
out to study the reactant transport and electrochemical reactions 
in a PEM fuel cell. Two-dimensional numerical analyses have 
been conducted along the flow direction by Hum and Li [10] and 
Chen et al. [11], and normal to the flow direction by Jeng et al. 
[12], among many others. In these works, the reactant transport 
was modeled including GDL on the assumption that the flow 
channel is straight and symmetrical to the boundary. Several 
three-dimensional simulations have been performed by Um and 
Wang [13], Hu et al. [14] and Hwang et al. [15] for the interdig- 
itated flow channel. They reported that the flow through GDL 
enhances mass transport and improve the performance due to the 
enhanced reaction area. Most of numerical modelings and sim- 
ulations, however, are based on symmetry considerations within 
the limited computational domains to reduce the computational 
costs [10-15]. Although such simplifications may be reason- 
able for simple reciprocating designs, e.g., parallel straight flow 
channels, it is not acceptable for the fuel cell with a long serpen- 
tine type flow channel. A typical serpentine flow channel has the 
cross section size in the order of 1 mm? or less and the length 
up to a few meters with many U-turns to fill the cell area on the 
bipolar plate as shown in Fig. 1. The pressure drop along the flow 
channel also results in a considerable pressure gradient across 
the GDL to the neighboring downstream channel due to much 
shorter distance. As a result, additional flow motion through 
GDL occurs, and is driven by the pressure gradient. This flow 
between the adjacent flow channels through the porous GDL has 
been referred to as flow cross over or cross leakage flow [16-25]. 
Convective flow motion through gas diffusion layer has been pre- 
dicted in 2D [16,25] and 3D [6,7,17,22] numerical simulations 
for a PEM fuel cell with conventional serpentine flow channels 
with various emphases. The effect of cross flow on cell per- 
formance was investigated through cell performance test using 
GDL with distinctively different gas permeability [19]. Dohle et 
al. [20] obtained flow homogeneity in GDL through flow visual- 
ization to find optimal permeability value for best performance. 

As aresult of this cross leakage flow, (i) the flow rate (or pres- 
sure) at any cross section of the flow channel cannot be treated as 
a known boundary condition, except at the inlet and outlet, since 
it is strongly coupled with this cross leakage flow; (ii) the effect 
of the cross flow is substantial on the total pressure drop in a 
cell because the cross flow brings about enlarging effective flow 
area. For example, present experimental results revealed that the 
amount of pressure drop in the presence of cross flow is up to 
80% smaller than that in the absence of cross flow for a typical 
PEM fuel cell of 100 cm? active cell area and 190 wm thick GDL 
with a single serpentine flow channel; (iii) the channel to channel 
symmetry, often used in CFD type simulation for the PEM fuel 
cell in order to reduce the computational domain, breaks down, 
as shown in ref. [16]. This symmetry breakdown requires the 
simulation of the entire cell instead of often employed in litera- 
ture and a much-reduced domain ranging from half anode flow 
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Fig. 1. Schematic of a polymer electrolyte fuel cell assembly. 


channel to half cathode flow channel with membrane electrode 
assembly sandwiched in-between; (iv) the most important effect 
of the phenomenon perhaps lies in the fact that the cross leakage 
flow through the GDL creates the convection flow through the 
porous electrode structure, resulting in effective water removal 
from the porous electrode structure. The cross leakage flow is 
significant [26] in actual PEM fuel cells, and needs to be included 
in any PEM fuel cell analysis and simulation, although it has 
been often neglected in previous PEM fuel cell studies. 
Therefore, the present study has been carried out to provide 
an experimental result documenting the degree of the cross flow 
effect on the reactant flow in a serpentine flow channel, and 
to provide a robust analysis on the cross leakage flow and its 
effects on flow distributions in a PEM fuel cell with a serpen- 
tine flow configuration. As pointed out early, three-dimensional 
numerical flow simulation of the entire flow channel, not a much- 
reduced domain, is necessary to capture the entire phenomenon. 
Hence, a three-dimensional numerical simulation has been per- 


formed covering the entire flow channel and GDL in a cell with 
the same physical dimensions and flow rates as the ones used in 
the experimental measurement. The effect of this cross leakage 
flow on the flow and pressure distributions has been analyzed for 
a wide range of flow conditions, as well as the permeability and 
thickness of the GDL. The numerical results are compared with 
the experimental measurements and good agreement is obtained. 


2. Experimental setup 


A PEM cell consists of a membrane electrolyte assembly 
(MEA) sandwiched between two bipolar plates as shown in 
Fig. 1. The gas diffusion layer (electrode backings), catalyst 
layer and polymer electrolyte membrane are referred to as MEA 
where current is produced. Fuel and oxidant are supplied to both 
sides of MEA through the flow channels on the bipolar plates 
producing electron in the anode catalyst layer and water in the 
cathode catalyst layer. The polymer electrolyte membrane con- 
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Fig. 2. Experimental setup. 
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Table 1 

Dimensions and conditions in the present PEM fuel cell simulation 

Parameter Value 

Bipolar plate 
Number of channel, Nen 1 
Number of U-turn, Mturn 50 
Channel width, w (mm) 1.0 
Land width, s (mm) 1.0 
Channel depth, d (mm) 1.1 
Channel length, / (m) 5.05 
Channel hydraulic diameter, dn (mm) 1.0476 
Cell active area, A, (cm?) 100.1 

Gas diffusion layer 
Outlet pressure, Pcou (atm) 1 
Permeability, K (m?) 107! to 1076 
Porosity, € 0.5 
Thickness, 5, (wm) 50-200 

Cell general 
Assembling torque (Ib in.~!) 50 
Flow rate velocity (m s7!) 3-30 
Inlet humidity (%) 100 
Temperature (°C) 20 


ducts the proton produced by hydrogen oxidation to the cathode. 
Fig. 2 presents a schematic for the measurement of pressure 
drop in a PEM fuel cell. Measurement was carried out for the 
cathode side since the cross flow is more significant owing to 
the higher flow rate in cathode stream, although measurements 
for the anode side can be measured with the same apparatus. 
A Fuel Cell Automated Test Station (FCATS S-800, Hydro- 
genics Inc.) was used for the present experiment. The FCATS 
system consists of a gas supply and conditioning sub-system, 
electronic load box, and a computer based control interface 
which performs data acquisition as well, a detailed descrip- 
tion of the FCATS is available elsewhere, e.g., [27]. The flow 
rate of the gas stream can be controlled within the +1.0% of 
the full-scale flow (16,000 sccm). Using the technique of steam 
injection, filtered de-ionized water is supplied to the gas streams 
to control the relative humidity. One hundred percent relative 
humidity air was employed in all of the experiments. Single cell 
test fixtures with active areas of 100.1 cm? were assembled for 
testing using MEAs manufactured by Ion Power Inc. All param- 
eters including flow conditions and dimensions are given in 
Table 1. 

The pressure drop from the channel inlet to the channel out- 
let was measured for two different conditions: in the presence 


Outlet 


a 


Inlet 


Fig. 3. A geometrical model of the flow channel and electrode backing layer. 


and absence of cross flow. For the latter case, the test cell was 
assembled with an impermeable film in place of MEA. Test 
cells were checked for gas leakage from the cell peripherals 
after the cell assembly and before the actual experimental mea- 
surements. Then before each measurement to be conducted, the 
cell assembly was pressurized up to 150kPa with the chan- 
nel exit closed. The cell was considered to have no leakage if 
the pressure was sustained for 1 min. Pressure data was stored 
with flow conditions after the system reached a steady-state 
which was typically within 1 min after a new flow rate was 
set. 


3. Numerical simulation 


The computational domain for the present numerical simu- 
lation includes a GDL and a serpentine flow channel, which is 
grooved on the bipolar plate, as shown in Figs. 3 and 4. Fully 
humidified air is considered as the reactant gas mixture in the 
cathode side and the following assumptions are made for the 
present simulations: 


(1) The gas mixture behaves like ideal gas, and is incompress- 
ible. 

(2) The GDL is considered to be isotropic porous media with 
constant porosity. 

(3) Flow is isothermal and laminar. 

(4) Water exists as a vapor only, hence only single-phase flow 
is considered. 


Polymer Electrolyte Membrane 
\ ZL Catalyst Layer 


d ~ Electrode Backing Layer 


Flow channel 


a 


Fig. 4. The details of computational domain and grid arrangement for numerical simulation. 
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Fig. 4 shows a typical three-dimensional structured grid for 
the flow channel and GDL used for the present numerical anal- 
ysis. In the flow channel, the flow field is obtained by solving 
the steady-state Navier-Stokes equations: 


V - (pu) =0 (1) 
V -(puu) = —V P + uV°u (2) 


where p is density, u the velocity, P the pressure and ju is the vis- 
cosity. The pressure drop in the GDL is proportional to velocity 
if the flow is laminar. This can be modeled by the addition of 
a viscous loss term according to Darcy’s law [16] such that the 
governing equations for the GDL becomes: 


V - (pu) =0 (3) 
Vy. (e=) = —eVP + uV°u — eu (4) 


where € is the porosity and K is the permeability of GDL. Since 
the pressure and Darcy’s terms are dominant over other terms 
Eq. (4) may be simplified in GDL, as often used in fuel cell 
analysis in literature: 


u 
VP= z u (5) 
where u is the superficial velocity in the porous medium. The 
constant pressure was given at the exit of the flow channel 
as boundary condition while the inlet boundary condition was 
assumed to have constant velocity. All other surfaces except the 
inlet and the outlet are assumed as no slip walls. 

The conservation of mass and momentum equations were dis- 
cretized using the finite volume method (second order upwind 
scheme) and solved by the commercial software FLUENT 6.0. 
It costs large amount of computational resource to simulate the 
entire flow channel with GDL since the flow channel and GDL 
has high aspect ratio of length to thickness. The length of the 
flow channel is 5.05 m with 50 U-turns on 10 cm x 10 cm of cell 
active area while the thickness of the channel and GDL are 1 mm 
and 50-200 um, respectively. At least more than 5 grid points 
were given for the channel width and the GDL thickness, which 
is equivalent to 200 um for flow channel and 10-40 um for the 
GDL, respectively. The number of the computational cells in the 
flow channel may seem a bit low at the first look, however, it is 
reasonably sufficient for the present purpose. Since this is vali- 
dated by doubling the number of computational cells around the 
bend region. Secondly, the significant cross flow shown in this 
study that injects the fluid into and out of the supposed secondary 
flow vortices zone around the channel bend areas weakens the 
secondary flow significantly. Similarly, Shimpalee et al. [22] 
used 6 computational cells for the flow channel in a related sim- 
ulation. A total of 1.5 million grid points have been generated 
for the entire computational domain consisting of the flow chan- 
nel and GDL. The Eqs. (1)-(4) were solved iteratively and the 
solution was considered to be convergent when the relative error 
in each field between two consecutive iterations was less than 
1074. Each case took about 7-10 h iteration time using the per- 
sonal computer equipped with dual CPUs and 4 Giga-byte of 
main memory. 


Table 2 
Bend head loss correlations [29] 
Configuration Correlations 
; y” 
Bending A Phend = ee 
2 k= Kioe T Ky 
ih- 0.21 
` kios = Tosa F Ja ps (90° bend, square duct) 
0/ah 


R 
kip = 0.0175 a (5 in degrees) 
h 


OR —(2/3) 
A= 32 (rea 4/ 2e) (Laminar) 
h 


Ro © 1.2dy in the present design. 


4. Pressure drop through engineering correlations 


The experimentally measured pressure drops from the chan- 
nel inlet to the channel outlet without GDL (impermeable walls) 
were also compared with the pressure drops determined by the 
Darcy—Weisbach’s correlation [28] in this study. The pressure 
drop caused by frictional loss along the flow channel can be 
estimated as from: 

L pv? 
AP = OF z (6) 
where Cf is the friction coefficient, L the length of the flow path, 
dp the hydraulic diameter of the flow path, p the fluid average 
density, and Vis the flow speed. The friction coefficient depend- 
ing on the flow condition is expressed in terms of Reynolds 
number [28]: 


a , Laminar: Re <2 x 10° 
Rea, 
Cr = (7) 
0.316 3 5 
Ta Turbulent : 4 x 10° < Re < 10 
J4 
Req, 


where Cı is 56.9 for the laminar flow in a square flow channel 
[28]. The Reynolds number is defined based on the hydraulic 
diameter of the flow path, 
_ PV 

H 


Re (8) 
where u represents average viscosity in the flow. Since the flow 
condition in flow channel is mostly laminar in the present study 
Eq. (6) can be reduced as: 


uL 
AP = OV (9) 
2dý 
The minor losses associated with the bending of the flow chan- 
nels are included in the estimation of the total pressure drop in 


the flow channel. Table 2 lists the minor loss correlations used 
in this study [29]. 
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Fig. 5. Flow structure in the first five flow channels and gas diffusion layer at several axial positions: (a) contour plot for the x-axial flow velocity, (b) two-dimensional 
velocity vectors (y and z component), (c) static pressure distribution. Vintet = 27 m s7!, ĝe =200 um, K=1.76 x 107!! m2 (K value is taken from [19]). 


5. Results and discussion 
5.1. Characteristics of cross flow 


The GDL, as porous electrodes for the anode and cathode 
side, is generally made of a porous carbon paper or carbon 
cloth, with a typical thickness in the range of 50-200 um in 
a PEM fuel cell. The porous structure of the electrodes enables 
the effective diffusion of the reactant gases to the catalyst layer 
in the MEA increasing the effective surface area available for 
the electrochemical reactions. The porous electrode also assists 
in water management by allowing the appropriate amount of 
water vapor to enter the MEA and humidify the polymer elec- 
trolyte. Although there exists a large amount of drag caused by 


porous structure the reactant transport can be driven not only by 
concentration gradient but also by pressure gradient in GDL. 
The flow structures for a typical serpentine flow channel and 
thin diffusion layer are shown in Fig. 5 in which (a), (b) and (c) 
represent the magnitude of axial velocity (x-component), two- 
dimensional velocity vectors (consisted of y- and z-component) 
and static pressure distribution at the three axial positions of 15, 
50 and 85 mm, respectively, for the first five channels from the 
channel inlet. Fig. 5(a) indicates that the axial velocity develops 
from the channel inlet and becomes developed gradually in the 
downstream direction with the maximum velocity near the cen- 
ter for the first channel. The axial velocity decreases drastically 
close to the GDL due to the substantial drag caused by the low 
permeability of the GDL. Itis interesting to observe that the max- 
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imum velocity as well as the cross sectional velocity distribution 
of the flow are skewed considerably in the second flow channel 
towards the downstream channel direction due to the cross flow 
from the upstream as well as to the downstream channels. It is 
seen that the degree of skew is increased for the third, fourth 
and fifth channels, or the successive downstream channels, 
because of the cumulative effect of the cross flow through the 
GDL. 

The flow structure in GDL is characterized by the two- 
dimensional velocity vector consisting of y- and z-component 
in Fig. 5(b). The cross flow velocity to the downstream chan- 
nel decreases as the flow approaches the U-turn (the end of the 
channel) where the pressure difference between the two adja- 
cent channels is the smallest. In contrast, the velocity from the 
upstream channel increases as the flow approaches the U-turn 
owing to increasing pressure gradient across the GDL. The static 
pressure distribution shown in Fig. 5(c) indicates that the pres- 
sure in the flow channel decreases downstream and that the 
lowest pressure gradient in the GDL occurs with the upstream 
channel before the U-turn and the highest pressure gradient with 
respect to the downstream channel occurs after the U-turn. It is 
also noticed that the pressure gradient in the channel area is 
much smaller than that in the GDL so that no pressure contour 
appears in Fig. 5(c) for the channel area. However, the pressure 
contour in the GDL is almost vertical in the area between the 
channels, indicating that the cross flow is from the upstream 
to the downstream channel with velocity almost parallel with 
the thickness of the GDL. However, the pressure contour in the 
GDL above the channel area is inclined, suggesting the veloc- 
ity is also inclined in this area, similar to the previous results 
[16]. 

Fig. 6 shows the average pressure in the GDL for various 
values of the GDL permeability and thickness which are the 
two major parameters affecting the amount and characteristics 
of the cross flow. As implied by Eq. (4), the permeability of 
the GDL has a direct and major impact on the cross flow, while 
the porosity will influence the cross flow indirectly through its 
relation with the permeability. The average pressure shown in 
Fig. 6 was obtained based on the two-dimensional pressure data 
set on the plane x = 50 mm for the GDL only and averaged for the 
same y positions. It is seen that pressure drops almost linearly 
in the y direction (in the cross channel direction) for all cases 
however the amount of pressure drop is considerably affected by 
the value of permeability in Fig. 6(a) and the thickness of GDL 
in Fig. 6(b). It is clear that the pressure drop in the cross channel 
direction, hence the cross flow between successive channels, 
is governed mainly by the permeability and thickness of the 
GDL. 

Fig. 7 provides four pathlines for the four fluid particles ini- 
tially located at four different vertical locations at the inlet flow 
channel cross section shown in Fig. 7(c). Fig. 7(a) shows the 
results for the GDL permeability of 107° m?, while Fig. 7(b) for 
the permeability of 10~!° m2. It is seen that the fluid particle a, 
being located on the top location, is the first being transported 
into the next flow channel by the cross flow, while the fluid 
particle d, being at the lowest position of the four investigated, 
travels the longest distance in the original flow channel before 


100 w 


80 "i 


AP (kPa) 


60 ere TS 


40 


(b) Y (mm) 


Fig. 6. Average pressure in the gas diffusion layer on the plane X = 50 mm in the 
y-direction: (a) 5g =200 wm, (b) K=107!° m2. The vertical coordinate AP rep- 
resents that the local average pressure is with respect to the pressure at the outlet. 


being transported to the next channel via the GDL by the cross 
flow. The left figures in Fig. 7(a and b) indicate that the fluid 
particles in the second, third, and further downstream channels 
are transported down the flow channel in the main direction of 
the channel flow before the cross flow takes them to the next 
downstream channels. However, the face view of the flow parti- 
cle trajectories shown in the right figures in Fig. 7(a and b) shows 
a striking difference. For the larger permeability GDL shown in 
Fig. 7(a), the fluid particles move from channel to channel by 
the strong cross flow without any noticeable convection effect 
of the main flow in the channels, except the first and last few 
channels near the inlet and exit. On the other hand, for the lower 
permeability results shown in Fig. 7(b), a much weaker cross 
flow results in the fluid particle trajectory significantly affected 
by the main flow in the flow channels. This highlights the com- 
plexities of the cross flow characteristics for different GDLs 
with various values of the permeability. Similar behavior can 
be observed as well for different thicknesses of the GDL as 
well. 
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Fig. 7. Pathlines illustrating the cross flows for four fluid particles starting from the four different start positions at the inlet of the flow channel: 5. =200 um, (a) 
K=10~° m’, (b) K=107!° m’, (c) the position of the four fluid particles at the channel inlet cross section. 


5.2. Quantity of cross flow; experimental measurement and 
numerical simulation 


The transport of reactant gas and water in GDL is the most 
critical aspect for interdigitated flow configurations [13-15] 
since no flow path is connected to the exit. Although the cross 
flow in PEM fuel cells with serpentine flow configurations have 
been reported before [16—18,22], its effect has been largely unac- 
counted for in most of designs and numerical modelings; for 
example, the amount of pressure drop was simply calculated 
using Eq. (6) including minor losses, assuming all the flow is in 
the channel. The amount of pressure drop in an actual fuel cell 
assembly has been measured for various flow conditions and 
compared with numerical results in Fig. 8. The PEM fuel cell 
for the present study has an active area of 100.1 cm? with one 
single serpentine flow channel (1.0mm width x 1.1 mm depth) 
and the thickness of the GDL was measured to be 190 um when 
the cell was assembled. Physical dimensions and flow condi- 


tions are listed in Table 1. The experimental results are given in 
symbols. It is seen that the pressure drop from the channel inlet 
to the channel outlet increases as the inlet Reynolds number is 
increased for both cases with the MEA sandwiched between the 
two flow plates (marked as “with GDL’) and with an imperme- 
able plate in between the two flow plates (marked as “without 
GDL’). It should be emphasized that the difference between the 
amount of pressure drop for two cases with and without the MEA 
is significantly large, as large as up to 80% of the pressure drop 
for the case without the MEA. 

In order to confirm that the significant difference in the pres- 
sure drop is indeed due to the cross flow phenomena described 
early, the results from both the three-dimensional numerical 
simulation and the engineering correlation, Eq. (6), including 
minor losses listed in Table 2, are also plotted in Fig. 8. It is 
observed that the results from the empirical correlation and the 
present numerical simulation are in reasonable agreement with 
the experimental measurements for the case without the MEA. 
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Fig. 8. Comparison of the amount of pressure drop vs. the Reynolds number at 
the channel inlet with and without gas diffusion layer (GDL): ôe = 190 um and 
K=1.76 x 107! m? used in the 3D simulation; experimental data from actually 
operating PEM fuel cell at the current density of 0.2 (+) and 0.4 A cm7? (x), 
stoichiometry (humidified air, cathode stream) of 3.0, operating temperature of 
80°C; all physical dimensions of the cell are kept the same as those given in 
Table 1. 


However, some difference exists at lower Reynolds numbers 
at the channel inlet. It is seen that for the low inlet Reynolds 
numbers the pressure drop estimated from the empirical cor- 
relation and numerical simulation agree with each other very 
well, although they both deviate from the experimental results 
with similar differences — this difference arises because in the 
experiments when the flow channel plates are assembled with 
the impermeable transparency replacing the MEA in the mid- 
dle, the transparency bulged into the flow channel area resulting 
in smaller flow channel sizes, hence higher pressure drop in the 
measured data. As the absolute pressure in the channel increases 
with the flow rate since the channel exit pressure is kept con- 
stant during the experiments, the transparency is being pushed 
back, the bulging effect is reduced and the prediction from both 
empirical and numerical results agree better with the experi- 
mental results as the flow rate or the inlet Reynolds number is 
increased. 

Fig. 8 also indicates that the experimental results for the 
case with the MEA are substantially lower compared to the 
case without the MEA, the present three-dimensional numerical 
simulation considering the cross flow phenomena again yields 
results that compare very well with the experimental results for 
the case with the MEA. Therefore, the substantial difference in 
the pressure drop with and without the MEA can be attributed to 
the effect of cross flow phenomena described in this study. Also 
in Fig. 8, two measurements are shown for the cell in operation at 
the current density of 0.25 and 0.5 Acm~?, corresponding to the 
cathode stoichiometry of 2. The water produced from reaction 
could be mostly in liquid form since the cathode inlet stream 
is saturated with water vapor. If the effect of the diffusion and 
electroosmotic drag is neglected, the resultant pressure drop for 
the running cell could be lower than that for the case with MEA 
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Fig. 9. Ratio of the volume flow rate for the cross flow through the gas dif- 
fusion layer to the inlet total flow rate Qinlet for various permeability values: 
Ointet = 3.08 x 1075 m? s™!, A-A’ is the surface only in the gas diffusion layer 
through which the cross flow rate Qanu is calculated. 


due to the consumption of the oxygen. However measured val- 
ues for the running cell are actually higher than the result for 
the case with the MEA in Fig. 8. This implies that liquid water 
might have existed inside the unit cell, that blocks the reactant 
flow, whether in the flow channel or in the GDL, resulting into 
higher pressure drop in the cell. 

From the preceding discussion, the effect of cross flow can 
be summarized as follows: The amount of pressure drop for an 
actual PEM fuel cell is significantly lower than the case without 
cross flow being considered. This lower pressure drop results in 
lower concentration gradient in cell with less pumping energy 
requirement. The reaction area can be enhanced by cross flow, 
which will bring about better performance, similarly to the cases 
studied for interdigitated flow configurations [13—15]. Further, 
the cross flow through the GDL could be an effective mechanism 
for the removal of water accumulated in the GDL. These factors 
together contribute to the better performance of PEM fuel cells 
using serpentine flow channel layouts, a fact that is well known 
for industrial PEM fuel cell developers and experimentalists [3]. 

The amount of cross flow is quantified as the volume flux 
through the first land area (between the first and second flow 
channel) and the results are presented as a relative amount to the 
total inlet flow rate in Fig. 9. As indicated previously, the cross 
flow is affected significantly by the permeability and thickness 
of the GDL. Permeability is a measure of the ability of a material 
(typically a porous medium) to transmit fluids through it. A typ- 
ical value of the permeability is within a range between 107!” 
and 10~!° m? for the GDL in a PEM fuel cell [6,15,30] while 
the thickness of the GDL may vary according to the designs 
and manufacturers. The amount of cross flow is considerable 
throughout the typical range of permeability and becomes dom- 
inant for the permeability greater than 10~!° m?. For instance, 
the amount of cross flow is more than 50% of the total inlet flow 
rate for the GDL thicker than 100 pm. For the MEA investigated 
experimentally in this study, corresponding to the results shown 
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Fig. 10. The total pressure drop from the flow channel inlet to the channel outlet 
for different thicknesses of the gas diffusion layer; K= 1.76 x 107}! m? taken 
from [19]. 


in Fig. 8, as much as 40% of the inlet flow does not follow the 
flow channel path, rather it crosses the land area through the 
GDL between the first and the second flow channel. Since the 
cathode GDL is most susceptible to liquid water flooding, this 
strong cross flow, coupled with a small width of the land area 
between the successive channels, would be beneficial for the 
removal of liquid water in the cathode GDL under the land area. 


5.3. Effects of permeability and thickness of GDL on the 
pressure drop 


It has been shown in the previous sections that the thick- 
ness and properties of GDL such as permeability are governing 
parameters influencing the cross flow in a PEM fuel cell with ser- 
pentine flow channels. The thickness of GDL can be measured 
after the cell is assembled, however porosity and permeability 
may not be known precisely for GDL as indicated previously. 
The permeability values for the GDL is quite different in liter- 
ature while the value for the porosity is typically in the range 
of 0.3-0.5 as noticed from refs. [12-15]. The strength of the 
flow resistance in the GDL, as shown in Eq. (4) or (5), is mainly 
influenced by the value of permeability and velocity. Therefore, 
simulation has been performed in a wide range of permeabil- 
ity and flow rate to investigate their effects on the amount of 
pressure drop in a cell. Fig. 10 shows the total pressure drop 
from the channel inlet to the channel outlet increases as the 
inlet Reynolds number is increased for all the GDL thicknesses 
investigated. However, the pressure drop for thicker GDLs is 
smaller due to a larger amount of cross flow through the GDL as 
shown previously in Fig. 9. The rate of increase for the pressure 
drop against the inlet flow rate becomes larger for thinner GDLs, 
attributable to the fact that the amount of the cross flow does not 
increase correspondingly with the inlet flow rate. 

Fig. 11 shows the total pressure drop for a wide range of 
permeability values for various thicknesses of GDL. It is seen 
that the pressure drop is larger for thinner GDL due to a weaker 
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Fig. 11. The total pressure drop from the flow channel inlet to the channel outlet 
for various permeability values of the gas diffusion layer: Vinlet = 27 m ah 


cross flow through such GDLs, and it decreases as the value of 
permeability is increased. However, the major decrement occurs 
within the range of permeability between 107!? and 1078 m?. 
The GDL behaves more like an impermeable wall for perme- 
ability lower than 107!? m? so that the pressure drop becomes 
almost identical for all different GDL thicknesses, and is equiv- 
alent to the pressure drop without GDL. On the other hand, the 
GDL provides additional flow path for larger values of perme- 
ability, so that beyond the permeability value of 1078 m? the 
pressure drop becomes almost independent of the permeability 
values, and is inversely proportional to the GDL thickness. 


6. Conclusions 


In this study the characteristics and effect of cross flow 
through the porous electrode structure between two adjacent 
flow channels have been investigated numerically and experi- 
mentally for a PEM fuel cell with a serpentine flow channel 
layout. Experimental measurements revealed that the total pres- 
sure drop from the channel inlet to the channel outlet in a PEM 
fuel cell is reduced substantially when compared with the case 
without the gas diffusion layer (that is, impermeable wall), and 
the amount of reduction can be as much as up to 80% for 
low flow conditions, although the percentage of reduction is 
decreased as the channel inlet Reynolds number is increased. 
Three-dimensional numerical simulation indicates that perme- 
ability and thickness of the gas diffusion layer are the two most 
important parameters influencing the cross flow and the resultant 
pressure drops. The total pressure drop is reduced for larger val- 
ues of the permeability and thickness of the gas diffusion layer, 
and significant reduction occurs within the permeability range 
of 107!? and 1078 m?. For the PEM fuel cell investigated in 
this study, a significant amount of the cross flow through the gas 
diffusion layer occurs, about 40% of the flow rate at the channel 
inlet. The present study suggests that cross flow through the gas 
diffusion layer between the adjacent flow channels is substantial, 
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and this strong convection current would help bring the reactant 
gas to the reaction sites and drive the liquid water out of the elec- 
trode structure for effective water management, hence partially 
responsible for the good performance of PEM fuel cells with ser- 
pentine flow channels. As a result, the effect of cross flow should 
be included for the analysis, simulation, design, operation, and 
performance optimization of practical PEM fuel cells. 


Acknowledgment 


This work was supported by AUTO21, the Network of Cen- 
ters of Excellence, Canada. 


References 


1] D.S. Watkins, K.W. Dircks, D.G. Epp. US Patent no. 5,108,849 (1992). 
2] D.S. Watkins, K.W. Dircks, D.G. Epp. US Patent no. 4,988,583 (1991). 
3] X. Li, Principles of Fuel Cells, Taylor & Francis, New York, 2005. 

4] X. Li, I. Sabir, Int. J. Hydrogen Energy. 30 (4) (2005) 359-371. 

5] V.A. Paganin, E.A. Ticcianelli, E.R. Gonzalez, J. Appl. Electrochem. 26 
(3) (1996) 297. 

6] J. Benziger, J. Nehlsen, D. Blackwell, T. Brennan, J. Itescu, J. Membrane 
Sci. 261 (2005) 98-106. 

7] G. Inoue, Y. Matsukuma, M. Minemoto, J. Power Sources 154 (2006) 8-17. 
8] H. Lee, J. Park, D. Kim, T. Lee, J. Power Sources 131 (2004) 200-206. 
9] L. Jordan, A. Shukla, T. Beehrsing, N. Avery, B. Muddle, M. Forsyth, J. 
Power Sources (2000) 250-254. 

[10] B. Hum, X. Li, J. Appl. Electrochem. 34 (2004) 205-215. 


[11] F Chen, Y. Wen, H. Chu, W. Yan, C. Soong, J. Power Sources 128 (2004) 
125-134. 

[12] K. Jeng, S. Lee, G. Tsai, C. Wang, J. Power Sources 138 (2004) 41-50. 

[13] S. Um, C. Wang, J. Power Sources 125 (2004) 40-51. 

[14] G. Hu, J. Fan, S. Chen, Y. Liu, K. Cen, J. Power Sources 136 (2004) 1-9. 

[15] J. Hwang, C. Chen, R. Savinell, C. Liu, J. Wainright, J. Appl. Electrochem. 
34 (2004) 217-224. 

[16] T. Kanezaki, X. Li, J.J. Baschuk, Cross leakage flow between adjacent flow 
channels in PEM fuel cells, J. Power Sources 162 (2006) 415-425. 

[17] S. Dutta, S. Shimpalee, J.W. Van Zee, Int. J. Heat Mass Transf. 44 (2001) 
2029-2042. 

[18] J.G. Pharoah, J. Power Sources 144 (2005) 77-82. 

[19] M.V. Williams, H.R. Kunz, J.M. Fenton, J. Electrochem. Soc. 151 (2004) 
A1617-A1627. 

[20] H. Dohle, R. Jung, N. Kimiaie, J. Mergel, M. Muller, J. Power Sources 124 
(2003) 371-384. 

[21] P.T. Nguyen, T. Berning, N. Djilali, J. Power Sources 130 (2004) 149-157. 

[22] S. Shimpalee, W.k. Lee, J.W. Van Zee, H. Naeri-Neshat, J. Power Sources 
156 (2006) 355-368. 

[23] P.H. Oosthuizen, L. Sun, K.B. Mcauley, Appl. Therm. Eng. 25 (2005) 
1083-1096. 

[24] Y. Wang, C.Y. Wang, J. Power Sources 147 (2005) 148-161. 

[25] W. Sun, B.A. Peppley, K. Karan, J. Power Sources 144 (2005) 42-53. 

[26] I. Sabir, Experimental investigation of proton exchange membrane fuel cell, 
MASc Thesis, University of Waterloo, March 2005. 

[27] A. Mughal, X. Li, Int. J. Environ. Stud. 63 (2006) 377-389. 

[28] F. White, Fluid Mechanics, third ed., McGraw Hill, 1994. 

[29] I. Idelchik, Handbook of Hydraulic Resistance, third ed., CRC Press, Boca 
Raton, 1994. 

[30] V. Gurau, H. Liu, S. Kakac, AIChE J. 44 (11) (1998) 2410-2422. 


